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We have calculated the electroweak radiative corrections at the O(a) level to 
the three channels of the process f\f\ZA —* and implemented them into the SANC 
system. Here A stands for the photon and f\ for a first generation fermion whose 
mass is neglected everywhere except in arguments of logarithmic functions. The 
symbol — > means that 4- momenta of all the external particles flow inwards. We 
present the complete analytical results for the covariant and helicity amplitudes for 
three cross channels: f\f\ — > Zj, Z — > /1/17 and fij — > f\Z . The one-loop scalar 
form factors of these channels are simply related by an appropriate permutation of 
their arguments s, t, u. To check the correctness of our results we first of all observe 
the independence of the scalar form factors on the gauge parameters and the validity 
of the Ward identity, i.e. external photon transversality, and, secondly, compare our 
numerical results with the other independent calculations available to us. 
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1 Introduction 



The group developing the network client-server system SANC (Support of Analytic and 
Numerical calculations for experiments at Colliders) actively continues to implement pro- 
cesses representing an interest for LHC and ILC physics. SANC is one of a few systems 
including Feynarts [1-3] and Grace-loop [4] in which calculations of elementary particle 
interactions were done at the one-loop precision level. A detailed description of version 
V.1.00 SANC was presented in Ref. [5]. The SANC client may be downloaded from two 
SANC servers Ref. [6]. 

In the recent papers [7, 8] we presented an ex- 
tension of the SANC Processes tree in the neu- 
tral current ffbb sector, comprising the version 
V.1.10. In this paper we realize its further ex- 
tension and include a calculation of the complete 
one-loop electroweak radiative corrections to the Z 
boson production channels ft ft — > Zy and fty — > 
fiZ, and to the Z boson decay Z — > fifty. This 
class of processes was already mentioned in section 
2.7 of Ref. [5] . For this reason, we do not change the 
number of SANC version, it is still V.1.10. The new 
processes are accessible from the ft ft — > ZA, Z — ► 
ftfiA and ftA — > f\Z nodes, which are placed 
in the Neutral Current sector of node 2f2b on 
the electroweak part (EW) of the Processes tree, 
see Fig. [TJ Each of these nodes contains standard 
modules of Scalar Form Factors (FF), Helicity Am- 
plitudes (HA), and bremsstrahlung (BR). 

The Zy production process is important for stud- 
ies of the anomalous trilinear Zyy and ZZy gauge 
boson couplings at the Fermilab Tevatron [9-11], 
LHC [12, 13] and at the Linear Collider [14, 15] in 
both the e + e~ and ey modes. The Standard Model 
(SM) of electroweak interactions predicts no trilin- 
ear gauge coupling of the Z boson to the photon at 
the tree level. Any deviation of the couplings from 
the expected values would indicate the existence of 
new physics beyond the SM. At the LHC, it is ex- 
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Figure 1: New processes in the ffbb 
sector. 

pected to observe hundreds of thousands of events 

of vector boson pair production. To match the precision of the LHC experiments, the 
vector boson pair production processes have to be considered beyond leading order [16]. 

Leptonic final states of the Z boson decays exhibit a very clear experimental signature 
and pave the way for precision tests of the SM beyond the leading order and possible 
detection of new physics. That is why it is necessary to fully control higher order EW 
corrections to the fermionic decays of the Z boson. 
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These processes were considered in the literature earlier mostly in connection with 
their sensitivity to anomalous triple gauge couplings, see for example papers [17-20]. 
To our knowledge, the QED and EW corrections to the Z boson production have been 
calculated previously only in papers [21-24]. 

All the processes under consideration can be treated as various cross channels of pro- 
cess /1/1Z7 — > 0, and hence one-loop corrected scalar form factors, derived for this 
process, can be used for its cross channels also, after an appropriate permutation of their 
arguments (s,t,u). This is not the case for helicity amplitudes, however. They are differ- 
ent for all three channels and must be calculated separately. 

The paper is organized as follows. In section [2] we demonstrate an analytic expression 
for the covariant amplitude at one-loop level in the annihilation channel. The helicity 
amplitudes for all three channels are given in section [31 In section H] we present numerical 
results computed by FORTRAN codes generated with software s2n and comparison with 
other independent calculations. Finally, summary remarks are given in section [51 



2 Covariant Amplitude 

Let us consider the process 

/i(pi, Ai) + /i(p 2 , A 2 ) + 7(p 3 , A 3 ) + Z(p A , A 4 ) -> 0, 



where the 4-momenta Pi (i = 1, 2, 3, 4) of all external particles flow inwards. Here, \{i = 
1,2,3,4) are the helicities of corresponding particles. Schematically this process is given 
in Fig. [2j where the black blob represents the sum of all tree and one-loop self energy, 
vertex and box type Feynman diagrams contributing to this process. The contributions of 
the counter term diagrams coming from the OMS renormalization procedure is assumed, 
as well. 




We found that next-to-leading order EW corrections to this process can be parame- 
terized in terms of 28 scalar form factors (FF) and corresponding basic matrix elements, 
14 vector and 14 axial ones. For the covariant amplitude (CA) we have: 
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where u(pi), u{p2) and are the bispinors and the mass of the external fermions, 
respectively; e3G°3) denotes the photon polarization vector and e^p^) is the Z boson 
polarization vector; the vector and axial gauge-boson-to-fermion couplings are denoted 
by Vf and a/, respectively; Ti a are the scalar FF of the vector and axial vector currents, 
respectively; F% and Str^ correspond to the lowest-order matrix elements. The usual 
Mandelstam invariants in Pauli metric (p 2 = —m 2 ) are defined as follows: 

(Pi+P2? = Q 2 = s, 



(p 2 +p 3 ) 2 = T 2 = -t, 

(P2+P4) 2 = U 2 = -U. 



(4) 



In Eq.(T5]) we keep the fermion mass in order to maintain photon transversality. Moreover 
in mass-containing denominators of Str° , the mass cannot be neglected because these 



denominators correspond to the propagators of fermions which emit external photons and 
thus would lead to mass singularities. 

The basic matrix elements, Str^,, are chosen to be explicitly transverse to the photonic 
4-momentum. That is, for all of them the following relations hold: 

Sti^fo), = 0. (5) 

We have checked that the FF Tl a are free of gauge parameters and of ultraviolet singu- 
larities (all calculations are done in the gauge). The analytical expressions of the FF 
are too cumbersome to be presented in this paper. They can be reproduced on-line with 
help of the SANC system. The CA for the processes we are interested in can be obtained 
from Eq.(j2J) exploiting crossing symmetry. This subject is covered in the next section. 

3 Helicity Amplitudes 

In this section we collect the analytical expressions of the helicity amplitudes (HA) for 
all three channels. Let us briefly recall the SANC strategy of observable (cross section, 
differential distributions) calculations. In a first step, SANC constructs the CA of the 
process, free of gauge parameters and of ultraviolet singularities, taking into account all 
lowest order and one-loop Feynman diagrams that contribute to the process. In the next 
step, HA are calculated analytically and converted into numerical code. Further, the cross 
section or the decay width of the process is formed as the incoherent sum of squares of 
all possible HA: 

da(dT)~ IWwjJ 2 **" (6) 

Ai A2...A n 

where squaring and summing is performed numerically. And finally, the Monte Carlo 
integrations over phase-space d$ n are performed using Vegas routine [25]. 

3.1 Annihilation channel f\f\ — > Z7 
To obtain the CA for the process 

fi(Pu Ai) + A(p 2 , A 2 ) -> 7(p 3 , A 3 ) + Z(p 4 , A 4 ), (7) 

where A»(i = 1,2,3,4) are the helicities of the external particles, we use the following 
substitutions of 4-momenta in Eq.(j2J): 

Pi -> Pi, 
P2 -> f>2, 
P3 -> -Ps, 
Pi -> -p 4 . 
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Figure 3: Schematic representation of one-loop Feynman diagrams for the annihilation 
channel. 



The set of non- vanishing HA for this process, which we denote as 'H\ 1 \ 2 \. i \ 4i1 read: 
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H ±T±T = =F-Z 4 (M Z ) sin$ 7 c+ 

o 

H ±TT ± = ±-Z A (M z ) sin$ 7 c_ 



_Z 1 (m f )Z 2 (m f ) 
2 
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with the following shorthand notation 
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(9) 



Here $ 7 is the center of mass system angle of the produced photon (angle between mo- 
menta f>2 and ps), t and u are the Mandelstam variables: 



t = m 2 f - Z 2 (m f ), 



u 
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(10) 



3.2 Decay channel Z — ► /1/17 

The CA of Z boson decay into fermion anti-fermion pairs and one real photon, 

Z(P2, A 2 ) -> 7(Pi, Ai) + /i(p 3 , A 3 ) + /i(p 4 , A 4 ), 
is obtained by interchanging of 4-momenta in Eq.flSJ) as follows: 
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Figure 4: Schematic representation of one-loop Feynman diagrams for the decay channel. 
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For the no n- vanishing HA, TC 
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where Tf and the coefficients k\ y2 are defined by Eqs. (J9]) and ( TTUT) with c± = 1 ± cost? 
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Here Z 2 (M Z ) = Mj — s and *&f is the angle between the vector p 3 and the direction 
defined by the photon momentum p x in the rest frame of compound (p3,p^)- The photon 
momentum, pi, is chosen to be direction of the z-axes in the (p3,Pa) rest frame. 



3.3 Z production channel e^f — > eZ 

And finally, in order to obtain the CA for the Z boson production channel 

7(pi, Ai) + e ± (p 2 , A 2 ) -> e ± (p 3 , A 3 ) + Z(p 4 , A 4 ) 
from Eq.(j2J), the 4- momenta permutations must be chosen as follows: 
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Figure 5: Schematic representation of one-loop Feynman diagrams for the Z boson pro- 
duction channel. 



The HA, H M A2A3A4, for this channel read 
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Here the coefficients &3,4,± are denned by 
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and t? e denotes the e ± scattering angle. The Mandelstam variables transform as follows: 
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4 Numerical results and comparison 



In this section we present the SANC predictions for various observables of all three processes 
under consideration. The tree level and single real photon emission contributions are com- 
pared with CompHEP, while one-loop electroweak and QED corrections for the production 
channel c-f — > eZ are checked against the Grace- loop package [4] and Ref. [24]. Note that 
all numerical results of this section are produced with the Standard SANC INPUT (section 
6.2.3 of Ref. [8]) if not stated otherwise. 

4.1 Annihilation channel fifi — ► Z7 

For this process we show in Table [1] a comparison between SANC and CompHEP results for 
the Born level cross sections and the cross sections of hard photon radiation. 





a, pb 


y/s, GeV 


100 


200 


500 


1000 


2000 


Born (SANC ) 


2482.0(1) 


86.230(1) 


11.652(1) 


2.9845(1) 


0.77816(1) 


Born (CompHEP) 


2482.0(1) 


86.230(1) 


11.651(1) 


2.9846(1) 


0.77817(1) 


Hard (SANC ) 
Hard (CompHEP) 


586.7(7) 
586.7(3) 


43.26(8) 
42.48(5) 


7.69(2) 
7.47(1) 


2.341(6) 
unstable 


0.717(2) 
unstable 



Table 1: Comparison of the Born and Hard cross sections of the e + e~ — > ^7(7) process 
(CompHEP input, > 1 GeV). The uncertainty of the last significant digit is given in 
brackets. 





a, pb 


y/s, GeV 


100 


200 


500 


1000 


2000 


Born (SANC ) 


1349.5(1) 


49.086(1) 


6.9785(1) 


1.8469(1) 


0.49555(1) 


Born (CompHEP) 


1349.4(1) 


49.086(1) 


6.9786(1) 


1.8469(1) 


0.49555(1) 


Hard (SANC ) 
Hard (CompHEP) 


173.82(3) 
173.82(3) 


14.138(3) 
14.083(21) 


2.7978(9) 
2.7627(21) 


0.9228(4) 
0.9045(11) 


0.3024(2) 
0.2936(4) 



Table 2: The same as Table [T]but for the process — > ^7(7). 

As can be seen from Tabled], we found very good agreement for the Born cross section. 
For the hard contribution we have perfect agreement at y/s =100 GeV, then a difference 
rapidly rising with energy, and eventually unstable CompHEP predictions for y/s at and 
above 1 TeV. As seen from Table [2] for the process — > £7(7), the hard contributions 
stay closer (though statistically incompatible) within a wider range of y/s pointing to the 
origin of the difference due to collinear singularities of the integrand. The stability against 
variation of uj discussed below gives us a great level of confidence in the SANC results. 
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In Tables I3HS1 we present the results of our calculations for the annihilation channels 
e + e~ — ► £7(7), uu — > ^7(7) and dd — > Zj(j), respectively, carried out with 10M 
statistics for the hard cross section for five energies and at each energy for two values 
of uj: uj = 10~ 5 \/s/2 (subscript 1) and 10 _6 v /i/2 (subscript 2); for a's in pb and for 

5 = ^l-loop/^Born _ j ^ % _ 



v/s, GeV 


200 


500 


1000 


2000 


5000 


o- Born , pb 


27.8548(1) 


3.37334(1) 


0.816485(2) 


0.202534(1) 


0.0323355(1) 


o-r iOOP , Pb 


43.36(4) 


5.216(9) 


1.239(4) 


0.299(1) 


0.0436(3) 


1— loop 1 

(7 2 p , pb 


43.38(5) 


5.211(10) 


1.235(4) 


0.298(2) 




5i, % 


55.7(2) 


54.6(3) 


51.9(4) 


47.4(6) 


34.9(8) 


5 2 , % 


55.7(2) 


54.5(3) 


51.3(5) 


46.9(8) 





Table 3: Comparison of the Born and one-loop cross sections of the annihilation channel 
e + e~ — > ^7(7) calculated with different values of the soft/hard separation parameter a); 
for details see the text. 



v/i, GeV 


200 


500 


1000 


2000 


5000 


a Born , pb 


4.7504(1) 


0.57540(0) 


0.13927(0) 


0.034548(0) 


0.005516(0) 


1— loop 1 

a l p , pb 


5.3399(8) 


0.6472(2) 


0.15367(6) 


0.036458(2) 


0.005203(6) 


1— loop 1 

a 2 p , pb 


5.3392(9) 


0.6470(2) 


0.17159(7) 


0.036458(2) 


0.005193(7) 


Si,% 


12.41(2) 


12.48(3) 


10.34(4) 


5.54(6) 


-5.67(11) 


62,% 


12.39(2) 


12.44(3) 


10.28(5) 


5.53(8) 


-5.84(12) 



Table 4: Comparison of the Born and one-loop cross sections of the annihilation channel 
uu — > £7(7) calculated with different values of the soft/hard separation parameter u. 



y/s, GeV 


200 


500 


1000 


2000 


5000 


o- Born , pb 


1.5230(0) 


0.18450(0) 


0.044658(0) 


0.011078(0) 


0.0017686(0) 


1— loop 1 

a l p , pb 


1.6033(1) 


0.18920(1) 


0.043823(4) 


0.009992(2) 


0.0012825(4) 


1— loop 1 

a 2 p , pb 


1.6033(1) 


0.18924(1) 


0.043825(5) 


0.009992(2) 


0.0012826(4) 


61,% 


5.274(4) 


2.549(6) 


-1.869(10) 


-9.807(14) 


-27.486(23) 


52, % 


5.275(4) 


2.570(7) 


-1.865(11) 


-9.804(16) 


-27.479(27) 



Table 5: Comparison of the Born and one-loop cross sections of the annihilation channel 
dd — > ^7(7) calculated with different values of the soft/hard separation parameter uj. 

The total 1-loop cross section ct 1_1oo p is the sum of the Born, virtual, soft and hard 
contributions: 

_1— loop ,^Born i „ virtual / \ \ i ^soft/x —\ i „hard/,— ,\ 

a 1 =0 + a {A) + cr [A,uJ) + cr [u). 
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Here cr virtual and cr soft depend on the regularizing parameter A which cancels in their sum. 
This cancellation was checked on the algebraic level. The contributions o" soft and <j hard 
depend on u), the soft/hard separation parameter. This dependence must cancel on the 
numerical level. To ascertain this cancellation we have done the calculation at each energy 
for two values of uj as shown above. Comparing the corresponding values of o- 1_loop and 
of 5 we can see that there is no change outside the statistical errors of the Monte Carlo 
integration. 

The following cuts were imposed: 

• CMS angular cuts for the Born, soft and virtual contributions where there is only 
one photon in the final state: i9 7i ^ G [1°, 179°] 

• CMS angular cuts on the Z boson and on the two photons and CMS energy cuts on 
the photons for the hard contribution: for the event to be accepted, dz and at least 
one of $ 7l or $ 72 must lie in the interval [1°, 179°], and both photons must have a 
CMS energy greater than uj. 

For all tables, the numbers in brackets give the statistical uncertainties of the last 
digit shown. 

4.2 Decay channel Z — ► /1/17 

In Table [6] we present the results of a comparison of the Born cross section and the cross 
section of hard photon bremsstrahlung of Z boson decay between SANC and CompHEP. We 
see that we have excellent agreement between these two programs. Differences are within 
statistical errors. 





T, GeV 


uj, GeV 


0.1 


1 


2 


5 


Born (SANC) 
Born (CompHEP) 


0.027730(1) 
0.027730(1) 


0.015779(1) 
0.015778(1) 


0.012269(1) 
0.012269(1) 


0.0078271(1) 
0.0078268(1) 


Hard (SANC) 
Hard (CompHEP) 


0.004393(2) 
0.004392(3) 


0.001358(1) 
0.001359(1) 


0.0007944(4) 
0.0007940(5) 


0.0002941(2) 
0.0002946(2) 



Table 6: Comparison of the Born and Hard widths of the Z — > [i + ^7(7) decay (Com- 
pHEP input, > uj for photon(s)). The uncertainty of the last significant digit is given 
in brackets. 

In Table [7] we show the differential decay rate dT/ds x 10 8 in GeV" 1 of the decay 
Z — > yU + /i"7(7), where \fs is the invariant mass of the pair calculated with two 

different values of the soft/hard separation parameter uj: 10~ 4 GeV (subscript 1) and 10~ 5 
GeV (subscript 2). The quantity 5 is given by 5 = (dV 1 - loop /ds - dT Bmn /ds)/dT Boia /ds. 

In Fig. [6] we show the differential decay widths dT Bmn /ds and dT 1 ~ loop /ds for the decay 
Z — > as functions of y/s = M M+M -. 

The Coulomb peak, which is due to photon exchange in the Feynman one-loop diagram 
with a 7Z7 three-boson vertex, is clearly seen. 
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y/s, GeV 


1. 


10. 


20. 


50. 


70. 


dT Boin /ds, GeV" 1 


7.918 


18.57 


22.63 


39.98 


89.11 


dY\- [oov /ds, GeV" 1 
dTl~ loop /ds, GeV" 1 


744.21(4) 
744.21(4) 


18.834(4) 
18.830(5) 


21.949(8) 
21.937(10) 


35.92(2) 
35.93(3) 


76.12(6) 
76.16(8) 


Si 

s 2 


92.992(5) 
92.992(5) 


0.0140(2) 
0.0137(3) 


-0.0300(3) 
-0.0305(4) 


-0.1014(6) 
-0.1014(8) 


-0.1458(7) 
-0.1452(9) 



Table 7: Comparison of the Born and one-loop differential widths of the decay channel 
Z — > /i + fj,~^('j) calculated with different values of the soft/hard separation parameter Co; 
for details see the text. 




Figure 6: Invariant mass distribution of the pair for the decay Z — > fi + fi 7(7). 

Both the Born (dashed line) and the 1-loop (dotted line) results are shown. 
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4.3 Z production channel cj — ► eZ 

As can be seen from Table [H we have again very good agreement between SANC and 
CompHEP predictions for the tree level and real photon emission cross sections of this 
process. 





a, pb 


y/s, GeV 


100 


200 


500 


1000 


2000 


Born (SANC ) 


82.266(1) 


23.716(1) 


5.5747(1) 


1.5343(1) 


0.40648(1) 


Born (CompHEP) 


82.265(1) 


23.716(1) 


5.5747(1) 


1.5343(1) 


0.40647(1) 


Hard (SANC ) 
Hard (CompHEP) 


4.012(1) 
4.014(0) 


3.689(2) 
3.688(1) 


1.368(1) 
1.364(1) 


0.4986(6) 
0.4973(6) 


0.1682(3) 
0.1678(3) 



Table 8: Comparison of the Born cross sections for the 'ye — > Ze reaction and of the 
Hard cross sections for the 7/x - — > Z^T'y reaction (CompHEP input, Ey > 1 GeV). 

In Table [9] we present the results of our calculations for the channel 7e~ — > Ze~{^{) 
carried out with 10M statistics for the hard cross section for five energies and at each 
energy for two values of a): Q — l§~^\/s/2 (subscript 1) and 10 _6 \/s/2 (subscript 2). 



s/s, GeV 


200 


500 


1000 


2000 


5000 


cx Born , pb 


8.3381(3) 


1.79168(0) 


0.46840(0) 


0.11842(0) 


0.019007(0) 


1— loop 1 

0"l , Pb 


8.7988(5) 


1.9591(2) 


0.52129(5) 


0.13171(1) 


0.02037(2) 


1— loop 1 

a 2 p , pb 


8.8002(9) 


1.9593(2) 


0.52131(6) 


0.13168(1) 


0.02037(3) 


61, % 


5.54(1) 


9.35(1) 


11.29(1) 


11.23(1) 


7.16(1) 


52, % 


5.54(1) 


9.36(1) 


11.30(1) 


11.20(1) 


7.15(2) 



Table 9: Comparison of the Born and 1-loop cross sections of channel 'je — > Ze (7) 
calculated with different values of the soft/hard separation parameter to. 

The notation of the various contributions, <7 Born etc., is as in the previous case. 

The cancellation of the A-dependent terms was again checked on the algebraic level. 
The cancellation of the u) dependence on the numerical level was tested as in the previous 
case. Comparing the corresponding values of cr 1_loop and of 5 we can see again that there 
is no change outside the statistical errors of the Monte Carlo integration. 

The following cuts were imposed: 

• CMS angular cuts for the Born cross section and for the contributions with Born-like 
kinematics: $ e , z G [1°, 179°] 

• CMS angular cuts on the Z boson and on the photon and a CMS energy cut on the 
electron for the hard contribution: for the event to be accepted, dz and d e must lie 
in the interval [1°, 179°], and the photon must have a CMS energy greater than Co. 

The numbers in brackets give the statistical uncertainties of the last digit shown. 
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In Table [TU] we show the comparison of the Born cross sections: the angular distri- 
butions da/d cos d e and the cross sections integrated over the given angular intervals, 
as well as the 1-loop EW corrections 5, produced by three programs: that of Ref. [24], 
Grace-loop [4] and SANC . 



y/s, Uev 







Ket. |z4j 


r< „„„„ i„„„ 
Kyi ace-loop 








20° 


-Born „u 
U , pD 


u.oyor 




u.oyouo 








X 07 
0, 70 


-5.9b 




-5.955b 






90° 


—Born ~^"U 

cr , pb 


U.6491 




r\ r? a ana 
0. 649116 








X 07 
0, 70 


O 

-O.OD 




O K ££0 

-o.55bz 


i nn 

1UU 




160° 


—Born 

a , pb 


n no o 




n no o o 








X 07 
0, 70 


1 n nn 




1 n nn cr. 
-lU.UUo 




20° 


< e < 160° 


—Born .^1^ 

a , pb 


i o n cr. i 
lo.Uol 




i o n cr i 
lo.Uol 








0, /O 


-y.U4 




n noon 
-y.Uooy 




1° 


< e < 179° 


—Born 

cr , pb 


on AO A 

33.484 




oo AO A 

33.484 








X 07 
0, 70 


i n 0^7 
-10.2 ( 




i n 070 
-10.2 ( o 






20° 


-.Born „U 
(7 , pD 


U.UZoyo 




n nosQS/i 

U.U/oyo4 








X 07 

0, /o 


on no 




on D7n 
-oU.U ( 9 






90° 


cr , pb 


n no c no 




n no c no c 








X 07 
0, 70 


-ZD. f4 




-zo. f 44 


500 




160° 


—.Born ~"U 

cr , pb 


n /i 1 
U.4obl 




n Afcczr\ r 7 
(J.4DDL)/ 








X 07 
0, 70 


oo nix 




oo n x /i 
-zo.Uo4 




20° 


< e < 160° 


—Born 

a , pb 


0.7051 


n Tn ci c 

0.70515 


n Tn cr i cr 

0.70515 








X 07 
0, 70 


-25.69 


o c c on 
-25. boy 


o c con 
-25.690 




1° 


< 9 < 179° 


—.Born 

a , pb 


1. ( (0 


1. (696 


1. (by ( 








X 07 

0, /o 


OO OT 


OO Q1 Q 


oo 01 o 
-ZZ.olo 






20° 


Born ^ 


001869 

V J > V 7 V 7 -L v 7 V 7 «.7 




0018688 

V 7 . V 7 V 7 _L V7 V./ v7 V 7 








5, % 


-41.57 




-41.575 






90° 


cr Born , pb 


0.002334 




0.0023340 








5, % 


-41.98 




-41.981 


2000 




160° 


CT Born , pb 


0.03094 




0.030942 








5, % 


-33.99 




-33.994 




20° 


< 9 < 160° 


CT Born , pb 


0.04620 


0.046201 


0.046201 








5, % 


-39.53 


-39.529 


-39.529 




1° 


< 9 < 179° 


Cr Born , pb 


0.1170 


0.1170 


0.11697 








6, % 


-30.84 


-30.845 


-30.845 



Table 10: Triple comparison of the Born cross section and of the correction 5 = 
^l-ioop / a Bom _ x for channel 1& - Ze~{i) (Ref. [24] input, E 7 = 0.025^ GeV). 

We have excellent agreement between these three results. Note that in this table the 
results taken from the literature were given there without statistical errors. The statistical 
errors of numbers obtained with SANC are in the digits beyond those shown. 
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5 Conclusions 



In this paper we describe the implementation of the complete one-loop EW calculations, 
including hard bremsstrahlung contributions, for the process fif±ZA — > into the SANC 
framework. The calculations were done using a combination of analytic and Monte Carlo 
integration methods which make it easy to calculate a variety of observables and to impose 
experimental cuts. We have presented analytical expressions for the covariant amplitudes 
of the process and the helicity amplitudes for three different cross channels: Z boson 
production /i/i — > and /17 — > f±Z, and for the decay Z — > /1/17. To be assured 
of the correctness of our analytical results, we observe the independence of the form 
factors on gauge parameters (all calculations were done in gauge), the validity of the 
Ward identity for the covariant amplitudes. We have compared our numerical results 
for these processes with other independent calculations. The Born level and the hard 
photon contrubutions of all three channels were checked against CompHEP package and we 
found a very good agreement except for the annihilation channel at high energies. For 
the channel 7e~ — > Ze~ ('-/), the comparison of the SANC EW NLO predictions with the 
results of Refs. [4,24] has shown an excellent agreement in a wide range of CMS energies 
and final electron scattering angles. 

The results presented lay a base for subsequent extensions of calculations in the anni- 
hilation channel appropriate to the process pp — > XZ^ at hadron colliders. 
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